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Abstract. We study the externally-driven motion of the domain walls (DWs) 
of the 7r/2 type in (in-the-plane ordered) nanostripes with the crystalline cubic 
anisotropy. Such DWs are much narrower than the transverse and vortex tt DWs 
in the soft-magnetic nanostripes while they propagate much faster, thus, enabling 
dense packing of magnetization domains and high speed processing of the many- 
domain states. The viscous current-driven motion of the DW with the velocity 
above lOOOm/s under the electric current of the density ~ lO^^A/m^ is predicted 
to take place in the nanostripes of the magnetite. Also, the viscous motion 
with the velocity above 700m/s can be driven by the magnetic field according 
to our solution to a ID analytical model and the micromagnetc simulations. Such 
huge velocities are achievable in the nanostripes of very small cross-sections (only 
lOOnm width and lOnm thickness). The fully stress driven propagation of the 
DW in the nanostripes of cubic magnetostrictive materials is predicted as well. 
The strength of the DW pinning to the stripe notches and the thermal stability 
of the magnetization during the current flow are addressed. 
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1. Introduction 

The domain walls (DWs) in the ferromagnetic nanowires (nanostripes) are the 
objects of a huge current interest because of concepts of miniaturization of the non¬ 
volatile memory and of the magnetics-based logic. Especially, the idea of the 3D 
digital recording known as the ” domain-wall racetrack memory” has taken enormous 
attention of the magnetics community mm- It is based on positioning a train of 
the magnetic DWs in the nanowire with the electric current. However, the racetrack 
concept has met difficulties. The most challenging of them is claimed to be a rapid 
increase of the temperature during the current-driven propagation of the DWs (a 
fast Joule heating) [21 H]. The minimum value of the spin transfer torque (STT), 
thus, the minimum intensity of the applied current and the minimum pace of the 
Jolue heating are limited by two factors. These are the DW pinning to natural 
and artificial (position-stabilizing) defects of the stripe and the need for reducing 
the operating time via increasing the DW velocity. Overcoming this difficulty has 
motivated the evolution of the DW racetrack architecture whose primary step was 
usage of the nanostripes with perpendicular magnetic anisotropy (PMA) instead of 
the longitudinally magnetized stripes [5]. However, since the PMA stripes are the 
uniaxial ferromagnets in principle, the field and current driven motions of the DW 
are turbulent whereas the record stability requires them to be viscous El IZ] ■ By the 
viscous motion we mean a uniform translational movement without the magnetization 
precession in the transverse plane (it takes place below the Walker breakdown). This 
requirement is satisfied in double-layer (or triple-layer) nanowires composed of a PMA 
ferromagnet (ferromagnets) and a non-magnetic metal. The stable (uniform) and fast 
current-driven propagation of the DW stack is possible utilizing the spin-Hall effect 
in the non-magnetic layer provided the global chirality of the DW system is stabilized 
by an interface-induced (Dzyaloshinskii-Moriya-like) anisotropy [5]. 

Simplifying the racetrack architecture is of importance with regard to the purpose 
of manufacturing the 3D memory device. The challenge is the stabilization of the DW 
structure in a simple nanowire (a single-layer nanostripe), thus, shifting the Walker- 
breakdown point upward the axis of the magnetic-field intensity. Our hypothesis 
is that the presence of two easy (or two hard) anisotropy axes could improve the 
stability. In particular, the racetrack created of a ferromagnetic material with a multi¬ 
axis crystalline anisotropy is an option to consider. Small width and high mobility of 
the DWs are desired and these advantages determine our preferences in terms of the 
material and sizes of the nanostripe. 

In a recent paper, analyzing the structure of the ferromagnetic DWs in the 
nanostripes made of cubic anisotropy materials, we (A.J. and P.G.) have predicted 
the tt/2 DW to stabilize in sufficiently-thin stripes (e.g. of Fe 304 ) provided the stripe 
axes coincide with the hard directions of the crystalline anisotropy [9]. Those DWs 
separate the domains magnetized along the easy axes in the stripe plane (diagonal to 
the stripe axes). Such an ordering requires the saturation magnetization to be small 
enough that the magnetostatic (shape-anisotropy) energy Ems = —fJ-o f m-Hd(m)dV 
is comparable to the exchange energy E^x = A fJ2i=x y z iVmipdH and it does not 
exceed the energy of the crystalline anisotropy Ean ~ Ki/M'^ J (writen 

up to the leading term). Here m denotes the magnetization, (M = \m\ and z axis 
is normal to the stripe plane), Hd denotes the dipole field, A and Ki < 0 are the 
exchange stiffness and cubic-anisotropy constant, respectively. In the static case, 
the cubic anisotropy is equivalent to the four-fold in-the-plane anisotropy due to the 
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Figure 1. A static 7r/2 head-to-head DW in the neighborhood of a head-to-tail 
7r/2 DW in: a tt; X r = 400 X lOnm^ cross-section nanostripe of the magnetite 
(a), and a static 7r/2 head-to-head DW in the neighborhood of many head-to-tail 
DWs in a tL! X r = 400 x 40nm^ cross-section nanostripe of the Terfenol-D (b). 
The colors and their intensity indicate the sign and value of the projection of the 
magnetization on the long axis of the stripe. 


planar alignment of the magnetization. The width of the 7r/2 DWs has appeared to 
be independent of the stripe width and much smaller than the width of the transverse 
or vortex DWs in soft-ferromagnetic nanostripes or in thicker nanostripes of the cubic 
materials as well. The small width of the 7r/2 DWs is expected to enable packing 
the bits (domains) in the racetrack with a density comparable to that in the PMA 
nanostripes. 

The aim of the present paper is to study the field-driven, current-driven, and 
stress-driven dynamics of the head-to-head (tail-to-tail) 7r/2 DWs in the nanostripes 
made of cubic-anisotropy ferromagnet with regard to the racetrack design. We 
estimate the velocities of the DW under different external drivings. The strength 
of pinning of the 7r/2 DW to a notch at one of the stripe edges is examined and the 
depinning current is determined. For the minimum current intensity established, the 
pace of the temperature increase due to the Joule heating is evaluated. Hence, we 
determine the main limitations of the 7r/2-DW based racetracks. 

The recent rapidly increasing interest in the methods of control of the 
magnetization state via the stress application (the straintronics) is motivated by 
the hope for reducing the energy consumption compared to the application of the 
magnetic field or the electric current. The strong crystalline anisotropy, in particular, 
the cubic anisotropy, is often accompanied by a significant magnetostriction [10] . This 
is shown below to offer a possibility of driving the head-to-head (tail-to-tail) 7r/2 DW 
propagation via the application of the mechanical stress. The study of the efhciency 
of such a driven motion completes the present analysis of the potential of the 7r/2 DW 
for the racetracks. 

In sec. II, we provide a background of the 7r/2 DWs in the cubic-anisotropy 
materials and consider the ways of driving the DW propagation. The field and current 
driven motions of the 7r/2 DWs in the nanostripe are investigated in subsec. Ilia, while 
the depinning from the notch is studied in subsec. Illb. Sec. IV is devoted to analyzing 
the stress-driven motion of the 7r/2 DW. Conclusions are collected in sec. V. 

2. Background 

In the previous work, we have performed the micromagnetic simulations of the 
DW formation in the nanostripes made of cubic-anisotropy materials and we have 
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developed an analytical description of the DW structures [5]. The density of the 
crystalline-anisotropy energy of the relevant cubic ferromagnets takes the form 
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see e.g. m- For a certain range of the material parameters (when the exchange 
lengths lex = \/‘2 A/plqM^, Ik = \J—AjKx are comparable), in very thin nanostripes 
[of the thickness of the order of the 7 r /2 DWs are stable solutions to 

the Landau-Lifshitz-Gilbert equation. In particular, they have been found with 
the material constants of the magnetite (Fe 304 ) for the nanostripe thicknesses of 
up to 20nm and with the parameters of the Terfenol-D (Tbo. 3 Dyo. 7 Fe 2 ) for the 
nanostripe thicknesses of up to 50nm, in the case of the crystallographic axes aligned 
parallel to the relevant stripe axes. The following material parameters of Fe 304 
have been included: M = 3 • I0®A/m, Ki = —l.I • lO^J/m^, K 2 = —0.3 • I0"^J/m^, 
A = 1.2 ■ 10“^^J/m, the Gilbert damping constant a = 0.014, and the gyromagnetic 
ratio 7 ~ • 2.21 • I0®m/As, [HI[121 [131 [M] . We remark that the effective saturation 

magnetization for the film of the magnetite is significantly smaller than for the bulk 
material (M = 4.8 • 10®A/m), which is due to the presence of anti-phase areas m- 
The mechanism of the anti-phase appearance is connected to an antiferromagnetic 
interaction at the anti-phase boundaries m- It is local and not expected to influence 
the energy of the single-ion anisotropy, thus, we use the anisotropy constants of the 
bulk. Let us mention that, contrary to previous reports on the attenuation of the 
anisotropy in the Fe 304 films, careful analysis suggests the four-fold in-the-plane 
anisotropy of the thin layers of the magnetite to be even stronger than in the thick 
films (the bulk limit) [16]. The relevant parameters of Tbo. 3 Dyo. 7 Fe 2 are following: 
M = 7.5 • 10®A/m, Ki = -8.7 • lO^J/m®, K 2 = 2.3 • lO^J/m®, A = 9.0 • lO-^^J/m, 
a = 0.1, [nUE]. 

The presence of the 7 r /2 DWs is a consequence of a large in-the-plane tilt of 
the direction of the domain magnetization from the long axis of the stripe, (see 
figure 1), that results from the strong crystalline anisotropy (in the case Ki < 0 ). 
The competition between the crystalline and shape anisotropies (the effect of the 
demagnetization field) results in a limitation on the length of the uniformly magnetized 
domains, thus, it leads to the creation of periodic structures of head-to-tail tt/2 DWs in 
the nanostripe (modulated superdomains). The superdomains of different projections 
of the magnetization onto the long axis of the stripe are separated by head-to-head 
or tail-to-tail 7 r /2 DWs, (similar structures of 7 r /2 DWs are observed in nm-thickness 
rings of Fe 304 , m)- Hence, unlike in the proposals of the tt-DW racetrack, we 
consider a record of bits which are written in patterned superdomains. The period of 
the patterns depends on the saturation magnetization (on the ratio plqM'^/K i^ where 
Po denotes the vacuum permeability) and on the thickness of the stripe since these 
two factors influence the demagnetization (dipole) field. 

Similar to case of the tt DWs in the stripes and wires, the uniform electric 
current through the cubic-anisotropy nanostripe induces STT in the DWs driving 
the propagation of all 7 r /2 DWs through the nanostripe in the same direction. Also, 
similar to the case of the tt DWs, the application of an external longitudinal magnetic 
field drives the neighboring head-to-head and tail-to-tail it 12 DWs to propagate in the 
opposite directions because the superdomains magnetized parallel (antiparallel) to the 
field must expand (shrink) [ 20 ] . 

Another way of driving the motion of the DW can arise when the neighboring 
ferromagnetic domains are not magnetized parallel to each other as it is in the 
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nanostripes under consideration. Any uniform in-the-plane stress applied in a different 
direction than the stripe axes is expected to induce a difference of the anisotropy 
energies of the neighboring domains due to the magnetostriction, thus, to drive 
expansion or shrinking of the domains In the case of a pair of one head- 

to-head and one tail-to-tail DWs magnetized in the same direction (of the opposite 
chiralities), the application of such a stress is expected to result in the mutual collision 
of the DWs and, provided the DW motion was viscous, in their annihilation. Hence, 
the stress application is a potential alternative to usage of the external magnetic field, 
in particular, with the purpose of manipulating the bit record via the annihilation of 
the DW pairs (the reversal of the superdomain magnetization). 

Because of the above, materials of strong (giant) magnetostriction are of our 
especial interest. Within the class of them, besides the magnetite and Terfenol-D, 
we have managed to simulate the 7r/2 DWs in thin nanostripes of the cobalt ferrite 
(CoFe 204 ). Remark that the nanotechnology of the Terfenol-D remains a challenge 
however. Unfortunately, simulating the domain wall creation in the nanostripes of 
the cubic-anisotropy and high-magnetostriction Fei_a;Gaa; (Galfenol-like) alloys (of 
the lOnm thickness at least), we have not managed to find the tt/2 DWs for any 
concentration x. This fact resembles quite high ratio of the exchange length of the 
crystalline anisotropy to the magnetostatic exchange length for Fei-ajGax, [9]. 

The density of the magnetoelastic energy of the cubic ferromagnet is of the form 


\ 1 2 2 , 2 2 , 2 2 \ 
Ai nn 777. /T -k 777, /T -k 777, /T 


n 


■me 



( 2 ) 


where a denotes the stress vector, and M = \m\ is the saturation magnetization 
m- Applied to the cubic media, one uses the name ’giant magnetostriction’ 
referring to the materials of large value of the volume magnetostriction parameter; 
A = 2/5Aioo + 3/5Aiii, where Aioo, Am determine the tetragonal and rhombohedral 
magnetostrictions respectively. The last formula holds for the polycrystalline 
ferromagnets |23] . When inducing the DW translation in a single-crystal film using 
the stress, the strong rhombohedral magnetostriction is needed. The energy of the 
tetragonal magnetostriction of the neighboring domains is independent on whether 
they are magnetized parallel or not provided the absolute values of the angles between 
the domain magnetization and the stripe axes are globally conserved. The energy of 
the rhombohedral magnetostriction makes a difference. Moreover, large value of the 
ratio Aiii/Aioo is an advantage (see Appendix B). Therefore, the materials of strong 
rhombohedral magnetostriction: Fe 304 , (Aioo = —19 • 10“® and Am = 81 • 10“®); 
CoFe 204 ; (Amo = -590 • 10“® and Am = 120-10-®) [H]; Tbo,3Dyo.7Fe2; (Amo = 
—90 • 10“® and Am = —1350 • 10“®) [17]; are of our especial interest. We notice 
that, in the contrast to them, the giant-magnetostriction Galfenol is of a near-zero 
rhombohedral magnetostriction. 

3. Field and current-driven domain-wall motion 
3.1. DW mobility 

In order to characterize the efficiency of the field-driven DW propagation and of the 
current-driven DW propagation, we determine the relevant mobilities of the head-to- 
head 7r/2 DW with micromagnetic simulations. In particular, the maximum velocity 
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Figure 2. Top: the dependence of the velocity of the 7r/2 head-to-head DW in 
the tL) X r = 110 X lOnm^ cross-section and w y. r = 300 x lOnm^ cross-section 
nanostripes of the magnetite on the longitudinal magnetic field (a), and on the 
current intensity (b). Bottom: dynamical transformations of the propagating 
head-to-head DW in the w y r = 300 x lOnm^ cross-section nanostripe of the 
magnetite under the the longitudinal field H = 1670A/m. The time step is 
At = 3.0ns, the reference frame is connected to the DW moving right the picture. 


of the viscous motion and the field intensity (current density) of the Walker breakdown 
are the characteristics which limit the applicability of DW-based devices. The velocity- 
field and velocity-current plots are shown in figure 2 a and figure 2 b, respectively, for 
the Fe 304 nanostripes of the widths r(;=110nm and w=300nm and of the thickness 
r=10nm. In the simulations (using the OOMMF package [24]), the cell size of the grid 
discretization is lOnm, while the stripe length is 30^m. In the plots, values of u(FI) 
and v(j) are average velocities over I5/rm propagation distance. It should be stressed 
that this propagation distance is sufficient to determine u(i/) in the considered range 
of the magnetic field and to determine v{j) in the regime of the viscous motion (below 
the critical current of the Walker breakdown). However, in the regime of turbulent 
motion (above the Walker breakdown), the propagation distance of 15/rm is insufficient 
to evaluate v{j) as an average over many cycles of DW transformations since the 
relevant propagation time is shorter than the period of the dynamical transformations 
of the DW. Therefore, the plotted in figure 2b v{j) curves do not resemble typical 
shape on the right-hand side of the local maximum (the Walker-breakdown peak) 
[25l [26]. On the other hand, for a reason discussed below (the fast Joule heating), 
there is no sense of simulating the DW motion on a wider time scale than the present 
one. 

For the purpose of analyzing the plotted data, a simple ID model of the field 
and current induced motion of the head-to-head (tail-to-tail) 7 r /2 DW is formulated 
in Appendix A. According to our analytical model, in the viscous-motion regime, the 
DW velocity depends on the intensity of the external field via the Walker-like formula 
v{H) = 'yp.oAH/a, where the DW width corresponds to the exchange length of the 
crystalline anisotropy A Ri y^—AfKi. Here, 7 denotes the gyromagnetic factor, ^0 
the vacuum permeability, a the Gilbert damping constant, A the exchange stiffness. 
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Basing on a 2D model of reference [5], the estimated in Appendix A width of the 
n/2 DW has been renormalized by the factor of two; A Ri 2y/—AjKi which leads to 
A r; 66nm for Fe 304 , and the DW mobility v{H)/H = 'ypoA/a = 1.04m^/As. This 
value coincides with the mobility determined from figure 2a for the tt/2 DW in the 
aspect-ratio w/t = 30 nanostripe. For the nanostripe of w/r = 11, we have established 
v{H)/H R! 0.55m^/As. Thus, the mobility of the 7r/2 DWs in the Fe 304 nanostripes 
is found to be at least two times higher than estimated and measured in reference 
[26] for the vortex DWs in the Py nanostripes of the aspect ratio w/t = 10 ^ 25; 
vpy{H)/H r; 0.2 ^ 0.3m^/As, (here and below, the parameters of the Py systems 
are indexed with the relevant superscript). Moreover, for the stripes of similar aspect 
ratio, the Walker breakdown fields in the magnetite structures are larger than in 
the Py structures. In the consequence, the maximum velocity of the field-driven 
viscous motion of the 7r/2 DW significantly exceeds that of the vortex DW in the Py 
nanostripes. 

Similar to the tt DW in the Permalloy, the head-to-head tt/2 DW undergoes a 
series of transformations into the vortex or/and antivortex DW above the Walker 
breakdown. There is characteristic peak in the field dependence of the DW velocity 
above the Walker breakdown that is well seen in v{H) plot for the nanostripe of 
the 300nm width (centered at the field value H tv 1.5 ■ lO^A/m) in figure 2a. It is 
accompanied by a transition between regimes of oscillatory and chaotic dynamical 
transformations of the DW following the nomenclature of reference m, where similar 
peak has been found for tt DWs in Py nanostripes. The transformation of the 7r/2 
DW into the antivortex DW is accompanied by the creation of two head-to-tail tt/2 
DWs on both sides of the head-to-head DW [see the magnetization snapshots for 
0 < t < 8 ■ 3ns in figure 2c], whereas, the transformation into the vortex DW and 
backward into the tt/2 DWs takes place during the passage of the head-to-head DW 
through a head-to-tail DW [see the magnetization snapshots for 8 • 3ns < t < 15 • 3ns 
in figure 2c]. 

Within the ID model of Appendix A, the velocity of the current-driven viscous 
propagation of the head-to-head (tail-to-tail) tt/2 DW is dependent of the current 
intensity following v{j) = fdijj/a. Here, rj is the constant of the STT strength; 
r] = PgpLB/2,eM, where, P, g, ps, and e denote the spin polarization, the Lande factor, 
the Bohr magneton, and the electron charge, respectively [26]. In the simulations of 
the Fe 304 nanostripes; P = 0.94, thus, 77 = 1.9 • 10“^°m^/C, [28]. The ratio fd/a 
determines the non-adiabaticity of the STT in the LLG equation [29]. Assuming the 
corresponding non-adiabaticity indicator in the LL equation to be equal to one, we 
take P/a = 2, (see [30]). Let us mention that, using the tt DWs in Py nanostripes 
(the polarization Ppy = 0.5, the saturation magnetization Mpy = 8.6 • 10®A/m), the 
relation jSpy/apy = 2 has been verified via measuring the DW velocity-to-current 
ratio vpy{j)/j = 0.73 • 10“^°m^/C, [3T]. 

We emphasize that current-related mobility of the head-to-head tt/2 DWs in 
the magnetite nanostripe; v{j)/j is higher for the magnetite nanostripe than for the 
head-to-head tt DWs in the Py nanostripe, (because of higher spin polarization and 
lower saturation magnetization). Also, the breakdown current of the tt/2 DW in 
the magnetite (the critical current of the transition between viscous and turbulent 
motion regimes) is at least twice the breakdown current for the vortex DW in the Py 
nanostripe [31] . Therefore, the velocity of the simulated current-driven propagation 
of the tt/2 DW is as high as llOOm/s for the current density of 3 • lO^^A/m^. It 
significantly exceeds the highest reported velocity of the DW in a sandwiched (triple- 
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Figure 3. The current-driven motion of the 7 r/ 2 -DW in the presence of a pinning 
site at the lower edge of the Fes O 4 stripe of the wxr = 110 x lOnm^ cross-section. 
The current densities are j = I ■ lO^^A/m^ (a), and j = 1.5 ■ lO^^A/m^ (b). The 
time step is At = 3.2ns. In b) the reference frame is moving relative to the pinning 
site. 


layer) nanostripe; 750m/s, that has been induced with similar value of the current 
density [32] ■ However, unlike the sandwiched nanostripe of reference ^2] (an artificial 
antiferromagnet), our ferromagnetic system enables the field-driven propagation of the 
DW as well! Moreover, following the present simulations, under currents exceeding the 
breakdown value, on a certain time scale that is shorter than the time of the dynamical 
transformation of the DW, (thus, the DW motion is thought of to be almost viscous), 
the velocity of the head-to-head 7r/2 DWs can be even higher than llOOm/s. 

The field-driven and current-driven motions of the head-to-head (tail-to-tail) tt/2 
DWs are different in another aspect than previously mentioned as well. Unlike the 
current, the longitudinal external field does not drive the head-to-tail tt/2 DWs to 
move. It is because the propagation of the head-to-tail DW is not accompanied by 
any change of the Zeeman energy in itself while the Zeeman energy changes during 
shrinking or expanding the superdomains (the propagation of the head-to-head and 
tail-to-tail DWs). In the consequence, the head-to-head 7r/2 DW propagates pushing 
one or two head-to-tail DWs. Any next head-to-tail DW that appears on the way 
of the complex of the three tt/2 DWs is annihilated during the collision, which is 
accompanied by a reduction of the complex to a pair of the DWs (one of the head-to- 
head type and one of the head-to-tail type). 

3.2. DW depinning and issue of thermal stability 

The issue to address is the strength of the pinning of the head-to-head (tail-to-tail) 
DW by stripe constrictions. The periodically distributed constrictions are considered 
to be pinning centers that stabilize the DW position in the racetrack counteracting 
mutual interactions between the DWs. We simulate the depinning of the DW from an 
artificially created (at one of the stripe edges) notch. The depth of the notch is 3/10 
of the stripe width while its shape is triangular. Remark that, in usual simulations 
of the current-driven motion of the DW, the depth of the notch is limited by the 
assumption of the uniform density of the current (the parameter of the STT is taken 
to be constant). That assumption could be contested when a variation of the stripe 
cross-section was large (the constriction of the stripe was very narrow). 

According to figure 3 and figure 4, the courses of the pinning/depinning processes 
are dependent on the orientation of the pinning-site triangle relative to the directions 
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Figure 4. The current-driven motion of the 7r/2-DW in the presence of a pinning 
site at the upper edge of the Fe 304 stripe of the w X t = 110 X lOnm^ cross- 
section. The current densities are j = 2.5 ■ lO^^A/m^ (a), and j = 5 ■ lO^^A/m^ 
(b). The time step is At = 3.2ns. 


of the magnetization of the closest domains. The strongest pinning is present 
when the the domains are magnetized parallel to the border of the notch, and this 
observation concerns the head-to-head (tail-to-tail) 7r/2 DWs as well as the head-to- 
tail 7r/2 DWs. The density of the depinning current for the head-to-head 7r/2 DW in 
w XT = llOnmxlOnm magnetite stripe appeared not to exceed 1.5- lO^^A/m^, which 
is several times higher value than achievable with PMA nanostripes while lower than 
the depinning current of the transverse DWs in the Py nanostripes of similar cross 
sections with similar notches [3]. Unfortunately, details of the notch shape are not 
standardized. A quite independent of those details indicator of the efficiency of the 
current-driven depinning is the ratio of the depinning current jp to the depinning field 

Hp. 

For the given depth of the notch, we have confirmed the theoretical relation 
jp = ap,Hp/l3r] with the simulations. Here p, = v{H)/H denotes the DW mobility. 
For instance, the ratio jp/Hp = 1.46 • 10^^Am“^/955Am“^ = 1.5nm“^ has been 
obtained from the simulations of w x t = llOnm x lOnm magnetite nanostripe with 
the geometry of figure 3 (the notch at the lower edge, the DW magnetized upward 
the stripe plane). It coincides with the theoretical value ap/fir] = 0.55m^C“^/2 • 1.9 • 
10“^*^m^C“^ = 1.5nm“^ Following reference [33], the maximum of the relevant ratio 
for the transverse DWs in the w x t = 300nm x lOnm Py nanostripe takes the value 
of jp^/Hp^ = 3 • 10^^Am“^/3500Am“^ = 0.86nm“^. Extrapolating this ratio to the 
nanostripe of lOOnm width (dividing it by three) leads to five times smaller value 
than for the head-to-head 7r/2 DW in the magnetite nanostripe {jp/Hp = 

Notice that the low efficiency of the current-induced depinning (high jp/Hp ratio) is 
not a disadvantage in itself. While the depinning current cannot be largely reduced 
with the magnetite nanostripe compared to that of the Py nanostripe, the depinning 
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Figure 5. The stress-driven motion of the head-to-head DW in Fe 304 nanostripe 
of the wXT = lOOx lOnm^ cross-section. The stress components are ax = 0.6GPa, 
ay = 2.4GPa, az = 0. The time step At = 0.5ns. 



field can be and it is found smaller than the Walker-breakdown value in the magnetite 
nanostripes. This enables the field-only driven viscous motion of the 7r/2 DW in the 
presence of the notches. In the case of the tt DWs, such a motion was impossible in 
the Py nanostripes nor with the PMA nanostripes. 

We note that, for a head-to-head DW magnetized upward the nanostripe, the ratio 
of the depinning current and the depinning field for the nanostripe with the notch at 
upper edge (hgure 4) is similar to that for the lower-edge patterned nanostripe of figure 

3. However, in the former case, the depinning current is only jp = 0.44 • lO^^Am”^. 
The strength of the mutual interactions of the head-to-head and tail-to-tail 7r/2 DWs 
requires a complex study (which is beyond the scope of the paper) in order to assess 
the minimal pinning field necessary. In particular, the mediating role of the structure 
of the head-to-tail DWs (whose period is scalable with stripe thickness) is to be 
determined. 

In order to determine the pace of the temperature increase due to the Joule 
heating, we apply the formula dT/dt = p/aCp, where, a denotes the electrical 
conductivity, C denotes the specific heat capacity, while p denotes the mass density. 
For Fe 304 nanostripe, in the presence of the current of the density j = lO^^A/m^, 
(cr = 1.0 • 10^1/Dm, C = 0.67J/gK, p = 5.0 • lO^kg/m^) [SH [35], one estimates 
dT /dt = 3.0- lO'^K/ns. The last value is significantly higher than for the Py nanostripe 
with the current of the same density [Ml [37]; dT/dt = 91K/ns, which is due to 
small electrical conductivity of the magnetite compared to Py. Notice that the Curie 
temperatures of Py and Fe 304 are similar, Tc ~ 850K. Any hope for slowing down 
the Joule heating of the magnetite nanostripe via material modifications (e.g. doping) 
in the direction of reduction of the electrical resistivity is weak since Fe 304 is the 
material of the lowest resistivity among all spinels. Hence, the limitation on the time 
length of the current pulses usable to drive the DW propagation in the magnetite 
nanostripes (~ lOps) is much stronger than for the Py stripes (~ Ins). However, the 
spatial distances between the neighboring head-to-head and tail-to-tail tt/2 DWs can 
be substantially reduced compared to the minimum spacing of the transverse/vortex 
DWs in the soft-magnetic nanostripes while the velocities of the DWs of the former 
type are an order of magnitude higher than those of the vortex DWs. Therefore, 
operating the Fe 304 structures with ultra-short current pulses can be as efficient as 
operating the Py structures with the ns pulses. 

4. Stress-driven domain-wall motion 

We verify whether the stress application is an efficient alternative to the field 
application when driving of the DW propagation. With this purpose, we have 
performed the simulations of the stress-driven motion of the DW for the magnetite 
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Figure 6. The mutual annihilation of two DWs (one head-to-head DW and 
one tail-to-tail DW) via the stress-driven collision in Fe 304 nanostripe of the 
w X r = 100 X lOnm^ cross-section, (ax = 0.5GPa, ay = 2.0GPa, az = 0, 
At = Ins). A transverse field of 30mT is applied in the direction (0,1,-1). 


stripe of the width w = lOOnm and thickness r = lOnm. The effects of the constant 
stress are easy to simulate with the MAGPAR package that we have utilized . The 
direction of the (in-the-plane) stress vector is taken with (Jxicfy = 1/4, Uz = 0, where 
X corresponds to the long axis of the stripe while z axis coincides with the out-of-plane 
direction. The stripe length is 1.5/im. 

With the snapshots in figure 5, it is seen that the dynamics of the DW consists 
of two stages. Upon switching the stress on, the ferromagnetic domains undergo an 
almost instant reorientation because the tetragonal magnetostriction creates novel 
easy directions. It forces the transition of the head-to-head (tail-to-tail) 7r/2 DW into 
the TT-like DW within Ins approximately. In fact, the stress-induced anisotropy due to 
the tetragonal magnetostriction is two-axial while both easy (in-the-plane) directions 
are slightly deviated from the long axis of the stripe. At the second stage of the 
simulations, the internal field due to the rhombohedral magnetostriction drives the 
resulting 7r-like DW to propagate along the wire. The propagation takes place under 
sufficiently high total stress a = \a\ > O.bvTlGPa. However, this lower bound on the 
motion-driving stress is not any fundamental value while the minimum stress changes 
with the simulation conditions (with the stripe length for instance). Its presence can 
be explained by changing the preference of the system in terms of the stabilization 
of the TT DW or another head-to-head DW (a n — eir DW, where e <C 1) under the 
stress. The later structure propagates while the former is static as mentioned in Sec. 
II. Turning on the stress in presence of a weak perpendicular magnetic field ensures 
the desired noncollinearity of the magnetization to be present. In particular, in the 
presence of a transverse in-the-plane field of 20mT, the simulated DW propagates 
under the (subcritical) stress of cr = 0.5-\/T7GPa. 

Notice that the stabilizing DWs are not in-the-plane magnetized while they are 
intermediate between the Bloch-like and Neel-like DWs. The stress-induced motion is 
almost viscous, however, minor oscillations of the DW structure are seen. 

The estimations based on the formulae of Appendix B, for the stress components 
(Tx = 0.6GPa, ay = 2.4GPa, az = 0 applied to the magnetite nanostripe (figure 5), lead 
to the DW width A m 26nm for cos((/7) = 1/10 (an almost Bloch DW). The velocity 
of the viscous propagation with dependence on the stress is given by formula (IB.31) . 
It leads, via fia = 270 • lO^m/sGPa, for the above values of the stress components, 
to the DW velocity 1600m/s that significantly overestimates the simulation result: 
210m/s (in the relevant figure 5, the stripe length is 1.5/J,m). The discrepancy is not 
a surprise, however, since the DW motion is not fully turbulent (slow) nor purely 
viscous (fast), as seen from figure 5. The intermediate (oscillatory) character of the 
motion is expected to be accompanied by a decrease of the average velocity compared 
to the estimated value. 
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In figure 6, we demonstrate (by means of the micromagnetic simulations) the 
stress-induced collision of the head-to-head and tail-to-tail 7r/2 DWs and the possibility 
of annihilating the DW pairs. Because of an instability of the direction of the DW 
magnetization, a weak transverse (to the long axis) field of 30mT is applied and it 
is deviated from the stripe plane by 7r/4 angle. Upon that stabilization, projections 
of the magnetizations of the colliding DWs are directed similar while the DW motion 
is viscous. Therefore, the expected result of the collision is mutual annihilation of 
the DWs following a general rule; the parallel (antiparallel) magnetized ferromagnetic 
DWs attract (repulse) [321 SO]- Our simulations confirm this expectation (figure 6). 

Addressing the way of application of the stress to the nanostripe, we mention 
the strain-mediated magnetoelectric effect in the multiferroic laminates. For instance, 
the magnetoelectricity of Fe 304 /PZN-PT and Fe 304 /PMN-PT is one of the strongest 
reported to date m- A rough estimation based on the data of reference [33] shows sub- 
GPa in-the-plane stress to be induced in a quite thick (sub-/im) structural magnetite 
layer upon the application of the electric field of 6kV/cm in the PZN-PT substrate. 
Moreover, increasing the stress in the magnetic layer via engineering the interface 
of the magnetoelectric composite is the subject of successive studies (e.g. m)- 

Therefore, we believe the desired values of the stress to be achievable with a simple 
voltage application. 

5. Conclusions 

We have analyzed the field, current, and stress driven propagation of the tt/2 DWs 
in the nanowires (nanostripes) created of cubic-anisotropy ferromagnetic materials, 
with relevance to the potential application of those textures to the DW racetrack 
memory. The tt/2 DWs are much narrower than the transverse DWs in the soft- 
magnetic nanowires of similar cross-sections which enables the recording with a bit 
density comparable to that of the PMA nanowires. In spite of the small width of 
the tt/2 DWs, the predicted velocities of the field-driven viscous motion of the tt/2 
DWs in the magnetite nanostripes are higher than achievable for tt DWs in the Py 
nanostripes. It is because of large values of the critical field of the Walker breakdown 
in the magnetite nanostructures. Let us recall that such a viscous field-driven motion 
is not achievable with the single-layer PMA nanowires. 

The current-driven propagation of the tt/2 DWs can be even faster than the fastest 
reported propagation of the chiral (Bloch) DWs in multi-layered nanowires with the 
PMA [33]. However, unlike for tt DW in the PMA nanowires, designing the racetrack 
for the tt/2 DWs, there is no need to involve other mechanisms of the propagation 
driving with current than the direct STT. The conservation of the global chirality of 
the system of many tt/2 DWs does not play any role in its positioning. 

The possibility of the stress-driven propagation of the DWs arises due to 
the anisotropy of the magnetostriction. We have demonstrated by means of the 
micromagnetic simulation the stress-induced collision of the DWs that results, under 
some stabilizing transverse field, in the annihilation of the DW pair. Hence, the stress 
application is found to be an alternative to using the longitudinal magnetic field when 
driving the domain magnetization reversal. 

The promising cubic-anisotropy material considered is the magnetite. First, 
because of a relatively low saturation magnetization (a weak shape anisotropy) and 
a strong cubic anisotropy, which enables the creation of the tt/2 DWs. Second, 
because of large ratio of the rhombohedral magnetostriction constant to the tetragonal 
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magnetostriction constant, which enables high-speed propagation of the DWs driven 
by the stress. Unfortunately, because of large electrical resistivity of the magnetite (the 
resistivity of the ferrites is two or three orders of magnitude larger than the resistivity 
of the Permalloy), the Joule heating remains a crucial problem in the racetrack design. 
However, small lengths of the domains and large velocities of the DWs enable as 
efficient current-driven positioning of the many-domain record as achievable with the 
Py-based racetracks, albeit using much shorter current pulses. 

Finally, we notice that the material parameters of the Terfenol-D seem to make 
it the most relevant material to manufacture the racetrack of the 7r/2 DWs. It is of a 
very strong cubic anisotropy (of a negative value of the Ki constant as desired) and 
of an extremely strong rhombohedral magnetostriction. A relatively small saturation 
magnetization (a weak shape anisotropy compared to the crystalline anisotropy) leads 
to the formation of very narrow tt/2 DWs in a wide range of thicknesses of the 
nanostripes (they are stable for up to 40nm-thick stripes). The electrical resistivity 
of the Terfenol-D is only 6.0 • 10“^Dcm while the specific heat capacity and the mass 
density are comparable to those of the magnetite and Py, thus, the Joule heating is slow 
[44] . A quite big (estimated) value of the Gilbert damping constant of the Terfenol-D 
is a potential disadvantage that limits the DW velocity however. On the other hand, 
the brittleness of the Terfenol-D makes it difficult to process to the nanoscale, (while 
the nanotechnology of the magnetite is advanced Hi]). 


Appendix A. Field (current) driven dynamics of 7r/2 DW in ID 


We treat a model of the tt/2 DW motion analytically solving the LLG equation 
dm a dm nj dm Bnj dm 

dt \ M dt dx 


qM dx 


(A.l) 


Tdca + Tdr. 


(A.2) 


in ID. Here, = —SH/Smi while 

Here, i = (1, 0,0), k = (0, 0,1), (the stripe is directed along the x axis while z axis is 
normal to the stripe plane), q denotes the gyromagnetic factor, A denotes the exchange 
stiffness. The longitudinal external field H is defined via B — {poH, 0, 0). The effective 
easy-axis (An) and easy-plane {K±) anisotropies are of the shape (magnetostatic) 
origin. The (cubic) crystalline anisotropy and the magnetoelastic coupling are included 
with the 'Hca-, 'Hme terms of the Hamiltonian, respectively. The constants of the normal 
and anomalous STT are denoted with rj and rj ■ fi, correspondingly, while j denotes 
the density of the electric current. 

Parametrizing the magnetization with angles; m = M[cos(d), sin(0) cos(</5), sin(0) sin(i^)], 
one describes a 7r/2 DW texture with 9 = 7r/4 + arctan(&/a), & = 1, a = exp{[a; — 
q{t)]/A{t)}. That ansatz for solving the equation (lA.ll) is a modified (via including the 
dependence of q and A on the time) static solution to the LLG equation [9|. Inserting 
the ansatz into (|A.ip . we notice that a/\/(a^ + P) = cos{9 — 7r/4), b/yj[a? -\- 6^) = 
sin(0 — 7r/4), thus, sin(d) = (a -b h)/^J2(a? + 6^), cos(0) = (a — h)/^2{a? + b'^). For 
qj = 0 (the in-the-plane ordering), and A = 0 (a sufficiently-long time upon switching 
the field or current on EilSD]), one arrives at the system of the equations 


a + m) 


ab 
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V ) {a? + 62)3/2 

- {llJ-oH + Lp) 

In the center of the DW [for x 


'yKw q 2 

= ( 7 ( 1 )], thus, for a = & = 


(AJ) 


Va 2 + 62 

1, these equations simplify to 


g + ??j 

2A 


+ a(/7 = 0, 


77Xoi? + tp 


aq + /3?7j 
2A 


(A.4) 


A rough estimate of the width of the DW using (|A.3|) leads to A = ^yA/{K\\ — Ki), 
where lATil » l/lOpoAl'^ > for the aspect ratio w/t ~ 10, (estimation details 
follow reference [26)1. 

The first of equations (IA.3I1 and the first of equations (IA.4I) become irrelevant 
to the description of the viscous DW motion since they originate from the dynamical 
equation of the out-of-plane magnetization component One removes them |25) . 
thus, finding q = {—Prjj + 'yAH)/a. 


Appendix B. Stress-driven dynamics of DW in ID 


We consider the head-to-head or tail-to-tail DW under an in-the-plane directed and 
sufficiently-strong external stress. According to the simulation of section IV, the 
domains which are magnetized in the (1,±1,0) directions in the absence of the stress 
rapidly remagnetize onto the long axis of the stripe upon switching the stress on. 
Thus, the head-to-tail DWs disappear while the head-to-head (tail-to-tail) 7r/2 DWs 
transform into tt DWs. The domain remagnetization is not complete however, while 
the magnetization direction is weakly deviated from the longitudinal one. That 
enables the creation of the longitudinal field inside the domains with the stress, thus, 
driving the DW motion can be explained with minimizing the energy via shrinking 
(developing) the domains magnetized almost antiparallel (parallel) to the field (similar 
to minimizing the Zeeman energy via the DW propagation). Additionally, an in-the- 
plane anisotropy is created with the stress, which influences the DW width. 

In order to model the stress-induced dynamics of the tt— e DW in the limit e — >■ 0+, 
we apply the standard ansatz 9 = 2arctan(6/a), & = 1, a = exp{[x — q{t)]/A{t)} 
that is used to describe the motion of the tt DW in the wires. Notice that sin(0) = 
2ab/+b'^), cos(d) = (a^ —&^)/(a^-|-6^). We consider the long-term dynamics, A = 0, 
in the absence of the external magnetic field and electric current. The substitution of 
the ansatz into the LLG system, for \ip\ <C 7r/2, leads to 


-I- = 0, 




27 [Ar||-I-Arj_sin^((/j)] 2'^Ki {of — b'^)'^ — Aafb'^ 27 A 


M 


37AiooCrCOs((73) 
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(72^-7?) 


b^-a^ 
a^ + b^ 


M {a? + b^Y A^M 

37Aiiicr7i72 (a^ — 6^)^ — 4a^5^cos^((/9) 


(B.l) 
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2ab{a? + 6 ^) 


= 0 , 


where 71 = 72 = o-y/cr, and tr = |cr|. In the center of the DW, a = 6 = I. 

Assuming the range of stresses K\\ + K±sm^{(p) <C —Ki — 3 AiooO'( 7 | — 7 i)cos((/?), the 
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A = 


2A 


-2K, 


3AiooO'(7| 


7 ?)cos((^) ■ 


(B.2) 


The last relation limits the range of (p angle to cos((/3) > 0, thus, the rotation of the 
DW magnetization, if takes place, is not stationary while oscillatory. We evaluate the 
velocity of the stress-driven (viscous) motion of the DW 

q = /ro-C0S^(v5)7i72cr, (B.3) 

_ 37A111A 
Ma ' 

Following (D, since tt DWs are completely out-of-the-plane magnetized, {p = ±n/2 
in the center of the tt DW), they cannot be driven by the stress, whereas, tt — e DWs 
which separate slightly noncollinear domains can, (ip ^ ±7r/2 in the center of such a 
DW). 
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